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HYDROMAGNETIC SHOCKS 
IN THE SOLAR WIND 

K. W. Ogilvie 

L. F. Burlaga 

Extraterrestrial Physics Branch 
Laboratory for Space Sciences 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 


ABSTRACT 

The Rankine-Hugoniot relations are applied to shock-like discontinuities 
measured by both magnetic field and plasma instruments on the satellite Ex- 
plorer 34 between May 30, 1967 and January 11, 1968. 

Shock normals were either determined from the magnetic field observa- 
tions, or from the times of occurrence of the discontinuity at Explorers 33, 34 
and 35. The Rankine-Hugoniot relations are obeyed to the accuracy of the ob- 
servations, and the values of shock velocities, density ratios, and Mach numbers 
indicate that at 1 AU the typical interplanetary shock is not strong, although all 
the events studied caused geomagnetic impulses. 
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HYDROMAGNETIC SHOCKS 


IN THE SOLAR WIND 


I. INTRODUCTION 

Shock-like discontinuities observed simultaneously in both the interplanetary 
magnetic field and plasma between May 30, 1967 and January 11, 1968 are ex- 
amined in this paper in order to determine whether they are consistent with the 
hydromagnetic theory for fast shocks. During the time when the satellite Ex- 
plorer 34 was outside the earth's bow shock and at a distance of more than 24 
earth radii from the earth, at least 7 shock-like discontinuities were observed. 
The selection criterion for such an event was the occurrence of simultaneous 
increases in the magnetic field intensity B, and the plasma density, bulk speed, 
and temperature. Since in a slow shock the tangential component of the magnetic 
field decreases across the discontinuity, these shock events , if they exist in the 
solar wind, would not be selected. All selected events were observed on terres- 
trial magnetometers as the sudden commencements of geomagnetic storms. 

The only previous observations of a propagating shock in the interplanetary 
medium, in which the discontinuity in the plasma properties and in the magnetic 
field were both measured, were made by Sonett et al. (1964), using detectors 
carried on Mariner 2. They reported observations of a discontinuous increase 
in the interplanetary plasma parameters and the magnetic field parameters 
which propagated with respect to the solar wind at a speed greater than the 
alfvdn speed. They showed that the velocity V x , density n x and temperature T 0 , 
(the subscripts 0 and 1 denote preshock and post-shock values respectively) pre- 
dicted from the Rankine-Hugoniot conditions for a fast, oblique, hydromagnetic 
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shock in an one-component, isotropic plasma were consistent with the observed 
values, but they found that the predicted temperature T x was 1.4 times higher 
than the observed temperature. Observations of events of a similar nature where 
either the plasma parameters or the magnetic field were studied, have been re- 
ported by Gosling et al. (1967a, b) Taylor (1968), Ness and Taylor (1968), and 
Van Allen and Ness (1967). Shock velocities have been deduced using conserva- 
tion relations by several of these authors. Taylor (1968) studied magnetic field 
observations taken by the satellite IMP-3 during 36 sudden commencements 
which took place in 1965, 1966 and 1967. He concluded that 26 of these events 
were caused by interplanetary shock waves, and after selecting eight of these 
with particularly well determined orientations, deduced that a typical shock front 
propagating from the sun to the earth has a radius of curvature somewhat less 
than 1 AU. This idea has also been advocated by Hirschberg (1968). 

H. EXPERIMENTAL 

The magnetic field observations used in this study were obtained by the tri- 
axial fluxgate magnetometer experiment of Fairfield and Ness. This instrument 
has a resolution of ±0.16y, and readings were taken every 2.56 seconds. 

The plasma instrument has been described by Ogilvie et al. (1968), so only 
a very short account will be given here. It records protons and helium nuclei 
separately, a spectrum of each species taking approximately one minute to ac- 
quire and successive spectra are separated in time by 3.04 minutes. The 
proton spectra alone have been used to deduce the fluid quantities used in this 
paper. This procedure might introduce errors of up to 20% into the calculation 
of the shock densities from (8), etc. However, the helium densities are not 
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- known on both sides of all the shocks, and this order of error exists in any case 
as a result of other uncertainties. 

The spin axis of the satellite is normal to the ecliptic plane to within 2.5°, 
and each of the fourteen energy per unit charge channels is sampled for 2.56 
seconds during which time the satellite makes a complete revolution. The ac- 
ceptance angle of the instrument is ±9° in a plane containing the spin axis, and 
2.5° in azimuth. 

ID. DESCRIPTION OF OBSERVATIONS AND 

SELECTION OF DISCONTINUITIES 

If it were known that the interplanetary plasma could be described as a one- 
component, hydromagnetic fluid with an isotropic temperature distribution, then 
the hydromagnetic theory implies that discontinuities across which density, bulk 
speed, temperature and magnetic field intensity all increase are necessarily 
fast shocks. In a multi -component, non-equilibrium plasma which contains about 
5% of helium ions, and is neutralized, as the solar wina is, by an electron com- 
ponent, it is still probable, but not certain, that such a discontinuity is a shock, 
assumed to be fast in all subsequent discussion. In any case, such a signature 
is a necessary condition for a shock. Thus, events for this study were selected 
by scanning ~3000 hours of interplanetary plasma and magnetic field data and 
identifying discontinuities across which n, V, T and B all increased by more 
than the respective errors in measurement. By a discontinuity, we mean that 
the plasma parameters changed in less than 3 minutes, (the instrument resolu- 
tion time) and the magnetic field parameters changed in less than 1 minute. Only 
events for which data v/ere available for several minutes before and after the 
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discontinuity were considered. Eight such clear discontinuities will be discussed 
in this paper. Other shock-like discontinuities were present (Burlaga and Ogilvie, 
1969) but the changes in the parameters were not large enough or sufficiently 
well known for the type of analysis which is presented below. The dates and 
times of the discontinuities are shown in the 1st two columns of Table I. The 
position of Explorer 34 and its distance from the earth at these times is also 
shown, labelled X, Y, Z and R in a coordinate system centered upon the earth, 
in which the positive X direction points to the sun, and the positive Z axis is 
northward normal to the plane of the ecliptic. 

Plasma Parameters at the Discontinuities 

The plasma instrument determines the differential proton flux in the 
streaming plasma at 14 energy values. From the non-zero members of these 
observations, we obtain 3 or 4 values of the distribution function , where v. , 

dv. J ’ 

) 

is the speed corresponding to the energy of step j . An approximation to the 
distribution function is obtained by piecewise fitting a maxwellian distribution to 
the measured dn/dv. (see Ogilvie et al., 1967). The fluid parameters are de- 
rived from the moments of this distribution function using the following equations , 

n = <v 0 > 


u = < v 1 >/<v 0 > 


T = k 


m <Vj> - <V,> J 


<V-> 


where 


V 
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These quantities were determined by averaging at 3 minute intervals for 
20-30 minutes on each side of each of the shocks in Table I. The pre-shock and 
post-shock values are presented in columns 2, 3 and 4 of Table I. 

In order to evaluate errors in density and temperature arising as a result 
of the fitting procedure, a computer program is used to simulate the process of 
measurement. Values of speed, density and temperature are deduced from the 
measured ’counts’, as discussed above. The ’counts’ which would have been ob- 
served by the detector in a plasma with a convected maxwellian velocity distri- 
bution having the same properties are calculated, and the fluid parameters re- 
determined from them. For the observations in Table I, it is found that the input 
and output velocities agree to ±5%, the densities to ±20% and temperatures to 
±50%, in the worst case. 

Magnetic Field Parameters at the Discontinuities 

Let B 0 be the magnetic field measured just before the shock arrived at the 
satellite and B x the field measured just after arrival. Since the magnetic field 
usually fluctuates appreciably near a shock, as discussed and illustrated for the 
case of the earth's bow shock by Fredericks and Coleman (1969) it is necessary 
to use some smoothing procedure to obtain B 0 and B x . The procedure used here 
is to average over six measurements (15 sec. in time) to compute B 0 and B x 
from the components. 

The magnitude, solar ecliptic latitude 8 rnd solar ecliptic longitude <p of 
the field vector obtained in this way for each of the events in Table I is shown in 
columns 6 and 7 of that Table. The last column of Table I shows a figure of 
merit which is defined by the equation 
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The first of these two factors characterizes the magnetic measurements. is 
the field change caused by the shock and ^I(SD) 2 j 1 ’ is the square root of the sum 
of the squares of the six standard deviations of the three field components on 
either side of the shock. These standard deviations are calculated for the 15 
second period over which the field values are averaged. Thus a small change in 
a varying field would give a low figure of merit. The second factor is obtained 
by forming the product of the density and the number of velocity intervals with 
non-zero flux. Thus a low density determined from a histogram with few bars 
would give a low figure of merit. 

IV. SHOCK NORMALS 

The normal to a shock surface, n, describes the orientation of that surface 
and also indicates the direction of propagation of the shock. As discussed in 
Section V, a knowledge of n is essential for a theoretical understanding of the 
changes which occur across a shock. This section discusses 2 methods for 
determining n , one using magnetic field measurements from one satellite and 
another using the shock arrival times at 3 satellites. 

Colburn and Sonett (1966) show that the shock normal vector may be deter- 
mined from the equation 

" = <B„ « 5,) x (B„ - x §,) « (B 0 - B,)| . (1) 

B Q and Bj are generally well enough known for (B 0 - B t ) to be well determined, 
but it is found that the magnetic field direction generally changes by only a 
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small angle due to the passage of the shock. The angle b In the equation 

|B, X B 0 j = Bj B q sin b 

varies between 0 and 25 degrees in the cases considered here, while the uncer- 
tainty in the magnetic field direction is 3 to 10 degrees, due to fluctuations near 
the shock. The accuracy of normal determination by this method was sufficient 
for three of the shocks, those on June 26, Sept. 13 and Sept. 19, and the results 
are shown in Table II and Figure 1. 

For the other events, where the angle b was insufficiently large compared 
to the field fluctuations, a second method, not previously used, will be described. 
It assumes only that the shock surface is planar over dimensions of order 50 R e . 
To illustrate this method, consider the discontinuity in the magnetic field on 
Jan. 11 observed by experiments on three spacecraft, Explorers 33, 34 and 35, 
which we shall refer to as 1, 2 and 3 respectively. 

We now refer to Table II , where we see that the spacecraft positions in the 
coordinate system used above were Rj = (58.8, 22.3, -28.4), R 2 = (17.9, -22.7, 
1.8) and R 3 = (-42, 46.4, 3.4). These three points define a plane P, and we 
transfer the origin in that plane to the position of satellite 1. Putting R 4 2 = 

Rj - R 2 , etc. we can define an angle u) by 

COS OJ - Rj 2 ■ ^13/1^12! ’ ’ 

and from Table n we see that R 12 = 93.9Re, R 13 = 108.4Re giving oj = 42.2°. 

The relative positions of the spacecraft in the plane P are illustrated in Figure 2. 
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If t 12 = (tj - t 2 ), etc. we seo from Table HI that for the Jan. 11 event 
t 12 = (8 ± 1) min., t 13 = (17.8 ± 1.8) min., and t 23 = (9.8 ± .8) min. The shock 
plane intersects the plane P in a line; let AB be a vector along that line, making 
an angle a with R 12 . From Figure 2 we can see that t 12 and t 13 depend on a , 
and that we define y to be the angle between the shock normal and the normal to 
the plane P. Thus U Sin y is the component of the shock velocity in P, and we 
can write. 


U sin y 


R 12 sin (a) 


t 


12 


R 13 sin (a; + a) 


t 


13 


Thus, 


t 12 R 13 

sin a - — ^ — sin (cu + a) , 

t 13 K 12 

and we find a = 30° ± 4° for Jan. 11. 

The vector AB is coplanar with R 12 and R 13 and can be expressed in terms 
of them. 


AB = -.74R 13 /R u + 1.42R U /R n 

so that AB = (.47, .84, -.22). We can now obtain the shock normal by forming 
the cross product of AB, a vector in the shock plane, with (Bj - B 0 ) 

n = ABx (Bj - B„). 

For the Jan. 11 shock this procedure gives a normal with a latitude angle of 2° 
and longitude of 153°. 
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Two other discontinuities, those on Aug. 29 and Nov. 29, were observed at 
the three satellites. The corresponding positions, time delays and the resulting 
shock normals and their latitudes and longitudes are given in Table II, and il- 
lustrated in Figure 1. 

V. RANKIN E -HUGO NIOT CONDITIONS 

Our aim here is to determine the extent to which the discontinuities in Table 
I can be described by the Rankine-Hugoniot conditions for a fast shock in a single 
component, magnetic plasma with an isotropic temperature. As has been pointed 
out by Krall and Tidman (1968), the distribution functions before and after a 
collisionless shock transition can be non-maxwellian, and there are many con- 
servation relations which could be tested in principle, beside the Rankine- 
Hugoniot relations. The nature, particularly the time scale, of the measurements 
prevents this at present, as they also prevent the study of shock structure. Ob- 
servations on the earth's bow shock indicate that a small flux of high energy 
particles, presumably produced by some acceleration process, move upstream 
from that discontinuity (Frank, 1968). These particles should be taken into ac- 
count in the detailed application of conservation laws to the bow shock, and 
similar effects might occur here. They have been neglected since they are un- 
likely to be the largest source of error. 

The plasma contains He ++ ions as well as protons; our observations of the 
helium are not precise enough to make comparisons with the theory for both 
components. We thus compare observations of the proton component alone with 
the single component equations given below. The average proportion of helium 
to hydrogen in the solar wind is 5% by number , while the ion temperature 
anisotropy is of order two. 
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The shock equations are as follows; 



(2) 


P V U^X 


B B 

II A 
477 


= 0 


(3) 


pv 2 + nkT + B u 


/877J - 


(4) 


. _ B, BB v 

y nkT " " M 


l A(v 2 + v 2 ) + — 1 — : — — + 1 „ 

' " 1 ' y - 1 P 477 p Air p Vj_ 


(5) 


The velocity components v x and v„ are defined with respect to a coordinate 
system centered at the shock and moving with the shock surface. The subscript 
1 refers to the direction normal to the shock surface, and the subscript II refers 
to the direction parallel to the shock surface and along the direction of the tan- 
gential component of B. Clearly these relations can be applied only to events for 
which n is known, since the equations involve the perpendicular and parallel 
components of v and B. 

Consider Equation (2). We measure the magnitude of the streaming veloci- 
ties | V Q | and | V x | before and after the shock in a frame which is stationary to 
the accuracy of the measurements with respect to the sun. Furthermore, these 
velocities are known to be radial, within a few degrees. Thus, the magnitude of 
the component of V along n, V-n s |v| cos x 0 , can be computed. Since the shock 
velocity is by definition Un, where U = |u| is the speed of the moving frame, 
v oi = I v 0 i cos X 0 ~ U* A similar expression may be written for v u , namely 
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v u = I coS ~ U. If it is assumed that the angle between V Q and Vj is 
small so that Xj = x 0 + e, where e<< 90°, then 

v lx = |VJ (cos X 0 + e sin X 0 ) " U 
v li = I Vj_ | cos X Q - U 


Because of the small (±9°) aperture of the plasma instrument, departures from 
radial flow of this order in a direction perpendicular to the ecliptic plane would 
be detected as gross decreases in particle density. Since such drops are not 
observed we can be confident that e is indeed small for the events in Table I. 
The expressions for V 0 and V lx substituted into 2 give 


U 


oypo v i - v o> 
(P \/ Pq - 1 ) 


cos X 0 . 


( 8 ) 


This is one necessary condition which must be satisfied if the discontinuities in 
Table I are shocks. Theoretical values of U, computed from this equation, are 
shown in Table III. No value of U is shown for the June 25 event, since an accu- 
rate normal direction cannot be obtained by either of the methods used for the 
other shocks. 


When U is known, a Mach number can be computed. We define the Mach 
number by the equation 


M 


2 

Ai 



This Mach number appears in the equations below. 

The theoretical shock speeds obtained from (8) can be compared with the 
observed shock speeds for the events of Jan. 11, Aug. 29 and Nov. 29. The 
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component of the shock speed in the plane P, illustrated in Figure 2 is given 
by the eouation 

R 13 

U sin y - ~ sin (co + y) . 

M3 


Now y is the angle between the shock normal n and n p , the normal to the 
plane P. This latter is (R 13 x R 12 )/ |R 13 x R 12 |. Writing n • n p = cos y, we get 
for Aug. 29, Nov. 29 and Jan. 11 respectively, y = 1°, y = 41° and y = 29°. The 
corresponding shock speeds are then (436*110) km sec -1 , (670±140) km sec" 1 , 
and (705±70) km sec -1 ; and these values are entirely reasonable by comparison 
with those obtained for these and other events by the application of mass con- 
servation, and illustrated in Table III. 


2.) Now consider Equation 3, which is a statement of conservation of 
momentum flux across a shock. Using (2), and the condition for continuity of 
tangential electric field in an infinitely conducting plasma, Wilkerson (1968) has 
obtained the following condition which must be satisfied if (3) is valid and if the 
discontinuities are shocks: 


n 


o 


n 


l 


B 


On 


B 


m 



(9) 


the ratio nj/n 0 was computed for all events for which a satisfactorily accurate 
shock normal could be obtained by either of the methods described above. When 
M a is greater than 2, the error in the calculated value of n x /n 0 is approximately 
20%; values of M A less than 2 lead to larger uncertainties of up to 40%. The ob- 
served and calculated densities agree within the stated errors, except for the 
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Nov. 29 event, for which the post-shock density was determined indirectly, 
using the relation between the change in the geomagnetic field and solar wind 
pressure discussed by Ogilvie, Burlaga and Wilkerson (1968). 

3.) Up to this point no assumptions have been made concerning the relation- 
ship between temperature and pressure in the plasma. Since the properties of 
the discontinuities examined above are consistent with magneto-hydrodynamic 
theory as far as we have gone, it is interesting to investigate Equation 4, which 
does not involve v n , and can therefore be solved immediately for the temperature 
ratio in terms of measured quantities. 



where 


r o - P 0 V <U /nkT O 31x1 ^0 ■ nkT 0 /B 0M /877 • 

This equation, due to Wilkerson (private communication) provides a test of 
Equation 4, assuming now a maxwellian velocity distribution on each side of the 
shock. The predicted temperature ratios are compared to the observed ratios 
in Table in. 

4.) Equation 5 cannot be verified because y is not known for the inter- 
planetary plasma. However, it can be used to determine an effective y which 
will guarantee that (5) is satisfied. 

Consider the relatively simple case when B ( %0. Equation 5 together with 
the other Rankine-Hugoniot conditions gives the well known relation (eg. 
Wilkerson, 1968), 
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where 


X r^/riQ - 1 + R(n 1 /n 0 - l) 3 
n i /n o ( X " n / n o) 


T/T 


o 


X 


7+1 

7 - 1 


Thus 


X = 


ni/n 0 (V T o - 1) 


( n i/n 0 ) 2 V T o " 


1 + 


|( ni /n 0 - l) 3 j 


The result for the discontinuity of Aug. 29 is 7 = 2. 3+1.0. The experimental 
errors are thus too large to distinguish between 7=5/3 and 7 = 2. 


SUMMARY 

We find as a result of examining the discontinuities detailed in Table I that, 
insofar as the difficulties associated with the small angular change in the mag- 
netic field could be overcome, their properties are consistent with being fast 
magneto-hydrodynamic shocks . 

The discrepancies shown in Table in, namely the predicted density ratio on 
Nov. 29 and the predicted temperature ratio on Jan. 11 are not too large to be 
the result of uncertainties in the determination of the shock normal. This work 
shows the determination of the shock normal to be by far the most difficult and 
critical part of the analysis of shocks , and emphasizes the advantages of multiple 
satellite observations. 


The only inconsistent event is that of Aug. 11, in which the only available 
method for determination of the shock normal is the application of the coplanarity 
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theorem. This should have given a well-determined result, but the value of 8 
obtained was 77°, inconsistent with other observed properties of the event. Ex- 
perimental difficulties associated with a flow direction at an appreciable angle 
to the ecliptic plane could have been responsible; it is not necessary to assume 
the event to be a shock whose properties would no* satisfy the Rankine-Hugoniot 
equations . 

The values of shock velocities, density ratios and Mach numbers indicate 
that at 1 AU the typical interplanetary shock is not strong, but the observations 
show that all these events caused geomagnetic storms. The Rankine-Hugoniot 
equations are obeyed to the accuracy of the data. 
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♦Post-shock density determined from geomagnetic fiei l change. 












Table III 


Calculated Properties of Shocks 


n i / n 0 ni/n 0 Tj/Tq 

Date & (n) <p (n) X Q U M A (B^) 

(obs) (calc) (obs) (calc) 


1967 Degrees km/ sec 


1 

June 25 

— 

— 

— 

— 

— 

2.3 ± .2 

B 

1.4 ± .4 

— 

June 26 

7 

357 

7.8 

482 

3.1 

1.9 ± .3 


1.2 ± .4 

1.7 

Aug 29 

8.5 

200 

21 

496 

15 

1.4 ± .2 

D 

1.9 ± .5 

1.8 

Sept 13 

-3.3 

4 

5.2 

416 

1.0 

1.5 ± .2 

1.0 

i 

1.4 ± .2 

1.6 

Sept 19 

-9.1 

335 

26.5 

497 

1 

2.1 ± .2 

1 

1.8 ± .3 

1.1 

Nov 29 

33.4 

163 

37 

394 


3.2 ± 1.0 


3.7 ± 1 

3.3 

1968 










Jan 11 

2.1 

151.4 

28.6 

524 

11.3 

2.6 ± .5 

1.5 

1.4 ± .5 

4.1 
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Figure 1. The distribution of the individual values of the elevations and azimuths of the six cal- 
culated shock normals. The dotted lines represent the multiple satellite observations and the 
sc 'id lines the directions obtained by application of the coplanarity theorem. Since all of the 
latter point away from the sun, the ambiguity in the former has been resolved in favor of the di- 
rection away from the sun. 
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